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ABSTRACT 



Analytical and numerical models are developed to predict 
the longshore sand transport distributions across the surf 
zone. The models, which also predict the root mean square 
wave height, H^, and the longshore current, V, are compared 
with field data acquired from Leadbetter Beach, Santa Barbara, 
California, during the intensive storm period in February, 1980. 
The breaker coefficient, B, and the bed shear stress coefficient, 
Cf, when equal to 1.1 and 0.005, were found to give the best 
agreement between the predicted H and V with the field 
measurements. The model for a plane sloping beach predicts 
the maximum sand transport at offshore distance X equal to 
0.8X^, where is the mean breakerline location, which agrees 
with Komar (1977a). The longshore sand transport formula, 
suggested by CERC (1977) , is used to calibrate the model which 
requires the empirical transport coefficient B g to be 0.18. 

The sand transport cross-shore distribution predicted by 
the model does not agree well with the corresponding trans- 
port inferred from the field measurements. Reasons offered 
to explain the differences are (1) that the model does not 
include the swash zone where the maximum transport is found 
and (2) that beach profile information beyond a distance of 
100 m offshore was incomplete. The results indicate the 
importance of including the swash zone to describe the effect 
of the long waves in predicting the sediment transport. 
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I. 



INTRODUCTION 



Water waves are one of the principal causes of shoreline 
changes. When waves break along the shore, they release their 
energy and momentum and give rise to a longshore current. 

The longshore current, along with the stirring action of the 
waves , is the primary mechanism for longshore sediment trans- 
port. The longshore sand transport rate is an essential fac- 
tor determining erosion or accretion along a coast. If the 
longshore sand transport rate can be accurately estimated, a 
quantitative picture of shoreline evolution can be evaluated, 
including changes of the shoreline due to marine structures. 

In the present study, analytical and numerical models are 
developed based on a longshore current model for random waves 
by Thornton and Guza (1983) and a sediment transport formula- 
tion by Thornton (1973) to predict the cross-shore sediment 
transport distribution and to compute the total volume of sand 
transport rate. The model is compared with the field data 
acquired from Leadbetter Beach, Santa Barbara, California 
( January-February , 1980) . 

The experiment at Santa Barbara was conducted as part of 
the Nearshore Sediment Transport Study (NSTS) . The objectives 
of NSTS are to develop improved techniques and engineering 
formulae to predict sediment transport on beaches having 
straight and parallel contours by utilizing field measurements 
of waves and currents at several nearshore sites. 
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During the Santa Barbara experiment, beach profiles were 
measured daily along with extensive measurements of waves and 
currents. A severe storm attacked the site of the experiment 
during 16-20 Feb., 1980 causing significant changes in the 
shoreline configuration and beach profiles. Despite the 
destruction of most of the instruments during the storm, 
a very substantial and valuable record of waves and currents 
was acquired which makes the Santa Barbara data very unique and 
significant for nearshore sand transport studies. 

The present work emphasizes the data acquired on 17 and 
18 February, 1980. A strong longshore current was observed 
due to storm waves approaching the beach at moderately large 
angles (6° at the breaker line) causing significant erosion 
in the beach profile. The daily beach profiles during the 
experiment ( January-February) show no well defined bars or 
troughs and essentially have a stable point 95 m offshore, 
which did not appear affected by erosion or accretion. Thus, 
it is hypothesized that the changes in the beach profile are 
due to sand transport in the longshore direction alone, and 
the on-offshore transport is negligible. Furthermore, the 
longshore sand transport distribution across the surf zone, 
predicted by the model, is assumed to have a distribution 
similar to the net longshore sand transport calculated from 
the profile changes. 
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II. LONGSHORE SEDIMENT TRANSPORT 



A. GENERAL 

When waves approach the shoreline at an oblique angle and 
break, they generate a longshore current parallel to the 
shoreline. This current in turn interacts with the waves to 
produce a longshore sand transport. The movement of sand 
particles is usually divided into two, not always distinct, 
modes of transport: bed load and suspended load transport. 

The bed load is the result of sand grains being rolled along 
the bottom by the shear of water moving above the sediment bed 
and is maintained by grain- to- grain contact. Bed load is 
usually considered to be confined to within 10 cm of the bottom. 
In the suspended load mode, the generally finer grains are 
transported by currents after they have been lifted from the 
bed by turbulent action of the wave-induced vertical motion. 

The total sand transport is the sum of the contribution by 
bed and suspended load transport. The rate at which the total 
volume of sand is transported parallel to the shoreline is 
termed "longshore sand transport: Q,". The two most success- 

ful and widely used approaches for estimating the longshore 
sand transport are an empirical wave power model and 
Bagnold' s energetics model. 

B. WAVE POWER MODEL 

In the wave power model, the rate of total longshore 
transport Q, is simply assumed proportional to the longshore 
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component of wave power (P^) at the breaker line to some power 
n (Watts, 1953a) . 




where and n are coefficients to be determined empirically. 
Assuming straight and parallel depth contours, the onshore energy 
flux in the direction of wave propagation per unit length of 
the beach is EC^ cos a, where E is the energy density of the 
waves, 0^ is the group velocity of the incoming waves, and 
a is the angle of the wave crests relative to the shoreline. 
Then, the longshore component of wave energy flux per unit 
length of the beach is given by: 

P„ = (EC ), cos a. sin a, (2.2) 

Ji g a d b 

where the subscript b refers to the breaker line. Applying 
linear theory. 



8 p 



2 

g H, C cos a. sin a, 
b g b 1 



(2.3) 



where p is the water density, g is the gravitational con- 
stant and H is the wave height. 

Watts (1953a) obtained the first field measurements by 
which the sand transport rate could be related to the local 
wave characteristics. Caldwell (1956) combined additional 
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data with those of Watts and found the coefficient (n) in 
(2.1) to be less than unity, and with a different proportion- 
ality coefficient (K^) than Watts. 

Inman and Bagnold (1963) interpreted field and laboratory 
data and obtained a linear relationship between Q £ and P^, 
i.e., n equal to unity, which takes into consideration the 
immersed weight of beach material. 



K 

(1-a) (P S ~P) 



(2.4) 



where p g is the density of sand transport, a is the porosity 
of beach sediment (taken as 0.40) and K is a dimensionless 
proportionality coefficient. 

Komar and Inman (1970) utilized sand tracers to measure 
longshore sand transport rates. Using this new data combined 
with data from Watts (1953a) , Caldwell (1956) and others, they 
obtained 0.77 for the proportionality coefficient K by relating 
the immersed weight transport rate I, to P, 



- °- 77 p * 



(2.5) 



3 3 

For quartz-density sand, p = 2.65 x 10 Kg/m , K = 0.77, 

s 

equation (2.4) becomes. 



Q, = 6.8 P, 



( 2 . 6 ) 
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3 

where Q in (2.6) is measured in units of m /day and in 
watts/meter. Komar and Inman use the root mean square wave 
height in (2.5) and (2.6) instead of the significant wave 
height, suggested by CERC (1977) , to evaluate the term. 

The relationships described in equations (2.5) and (2.6) 
are almost purely empirical with no real consideration of 
the sand transport mechanism. In addition, the coefficient 
0.77 in equation (2.5) has a high degree of uncertainty 

Longuet-Higgins (1972), in attempting to rationalize (2.2), 
points out that since the energy flux (EC^) is a vector rather 
than a second order tensor, the longshore component of wave 
power (P ) can be written equivalently as the product of two 
physically meaningful quantities: 

C 

P = (E cos a sin a]C = S C (2.7) 

36 ^ xy 



where S is the longshore component of the radiation stress 

and C is the phase speed of the waves. The radiation stress 

(S ) is shown by Longuet-Higgins (1972) , among others, to be 

the principal driving force for longshore currents, and thus 

would be expected to be important in describing longshore 

sediment transport. Vitale (1981) compared the radiation 

stress (S ) and the wave power (P 0 ) as predictors of the 

longshore sediment transport rate by conducting three-dimensional 

movable-bed laboratory tests and he found no major difference 

in choosing S over P to predict the longshore sediment 
xy Jo 

transport . 
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c. BAGNOLD'S ENERGETICS NO DEL 

Bagnold (1963) relates the rate of immersed weight trans- 
port (1^) to the work done by waves and current. The mechanics 
of sand particle movement are described as a back-and-f orth 
motion under the action of the wave orbital motion with 
essentially no net transport. Wave energy is expended in 
supporting and suspending the moving sand above the bed. 

Once the sediments are suspended, the presence of a unidirec- 
tional current superimposed on the to-and-fro motion can then 
produce a net drift of sediment. Bagnold (1963) derived the 
relationship : 



K' 



w 




u 

o 



( 2 . 8 ) 



where i n is sediment transport per unit width in the direc- 
tion 0 determined by the unidirectional current u^ , o is the 
available wave power supporting the sediments above the bot- 
tom, u q is the orbital velocity of wave motion, and K' is a 
dimensionless coefficient. 

Inman and Bagnold (1963) specified the various parameters 
of waves at the breaker line and obtained: 



v, 

I, = K' (EC ) , cos a. 
x, g b u b 

3 m 



(2.9) 



where u is the maximum horizontal orbital velocity of the 
m 2 

waves. Equation (2.9) can be viewed as a general relationship 
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because it does not specify the cause of the longshore cur- 
rent v^; it can be either tidal, currents of a cell circula- 
tion, wind generated, or due to oblique wave approach. Komar 
and Inman (1970) utilized their littoral drift measurements 
to test equation (2.9) and required K' to be equal to 0.28 
to agree with their data. In terms of and for quartz 
density sand, equation (2.9) is given by 

v, 

Q 0 = 2.5 (EC ). — (2.10) 

g b u 
^ m 

3 

where is in m /day and the wave and current parameters 
are in mKs units (watts/meter) (Komar, 1983). Wang and Chang 
(1979) found good agreement between their measurements along 
the bayshore of a barrier island in the Gulf of Mexico, and 
equation (2.9) for a coefficient K' = 0.18, where the long- 
shore currents generally are due to local winds and tides. 

Bowen (1981) applied Bagnold's equation to the problem of 
on-offshore sediment transport on beaches and then studied 
the special case of a normal incident wave and steady on- 
offshore currents (no longshore currents) . The waves are 
described using Stokes second order wave theory and the mean 
on-offshore current using Longuet-Higgin ' s (1953) bottom 
streaming solutions. The model predicts an equilibrium 
beach profile as a function of the incident wave character- 
istics. Bowen's results support observations that steep 
beaches are generally coarse grained, and short period waves 
cause a beach to erode. 
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Figure (2.1). Schematic of the Bagnold model for sand 

transport wherein the orbital velocity 
uq due to the waves place^ the sand in 
motion, and the current u^ provides a 
net transport of sand, I ^ . 
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Bailard and Inman (1981) independently derived a very 
similar model to Bowen's (1981) to describe the instantaneous 
transport in the nearshore zone, including both on-offshore 
movements as well as longshore movement. An important part 
of the study is the effect of local bottom slope on sand 
transport rate, as this permits the analysis of the formation 
of the longshore bars and other inshore topography. Bailard 
(1982) , extended the model to predict the total load trans- 
port for time-varying flow over an arbitrarily sloping planar 
bed. The model predicts the local equilibrium beach slope 
as well as the local near bottom sediment transport rate as 
a function of the near bottom velocity. Surf zone conditions 
are described as steady longshore and on-offshore currents 
in addition to a wave induced current having a local wave 
angle . 

A number of other models have been proposed. For instance, 
Madsen and Grant (1976) adapt the Einstein-Brown sediment 
transport equations to time-varying sand movements under com- 
bined waves and unidirectional currents. The bottom stress 
(t) under the combined oscillatory motions of the waves and 
any superimposed currents, and the resulting sand transport, 
vary with time. The model is relatively complicated and re- 
quires parameters for which no data are available. The 
success of the wave power formulation and Bagnold's energetics 
model is due to their relative simplicity requiring only one 
or two parameters specified from the data. 
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In general, the predictions of the littoral drift have a 
high degree of uncertainty. This partly results from our 
basic inability to make accurate measurements of sand trans- 
port on beaches. However, the two approaches of equations 
(2.5) and (2.9) can be used to give at least a qualitative 
description of littoral drift. 

D. DISTRIBUTION OF LONGSHORE TRANSPORT ACROSS THE SURF ZONE 
Bagnold's (1963) approach has been also utilized to pre- 
dict the distribution of longshore sand transport. Thornton 
(1973) analyzed sand transport inside and outside the surf 

zone where the volume transport rate per unit width (q ) is 

s 

given by: 



B 



g ( 1 - 7T-) 



1/2 3(EC g ) 

u 1 3x 
m 



( 2 . 11 ) 



where B^ is an empirical dimensionless coefficient. Thornton 

(1973) obtained field data at Fernandina Beach, Florida, to 

evaluate B in equation (2.11) . The transport rate (q ) was 
s s 

measured using a series of bed load traps placed in a line 
across the nearshore. The results indicate that the maximum 
longshore transport is in the breaker zone with B s = 0.08 
and stronger transport occurs over the bars rather than in 
the troughs (Fig. 2.2). 

Komar ( 19 76b , 19 77a) similarly utilized Bagnold's (1963) 
model to analyze the distribution of the littoral drift in 
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(2.2) . The distribution of the longshore send 
transport measured by Thornton (1973) 
at Fernandina Beach, Florida, by bed 
load traps in the positions as shown. 
The dashed lines are based on equation 
2*11 with B = 0.08, [from Thornton 
(1973)] s 
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terms of (I^). The local immersed-weight of sand transport 
rate per unit of width inside the surf zone (i x ) is given 
by : 



i(x) 




(0.5 C f ) 




h (x) v ( x) 



( 2 . 12 ) 



where C^. is the drag coefficient of oscillatory wave motions 
and h and v are the local water depth and longshore current, respec- 
tively; both are functions of the distance x offshore from the shore- 
line. K., is a dimensionless proportionality factor deter- 
mined by calibrating the model equation (2.12) to yield 1^ 
calculated from equation (2.5) where the total immersed- 
weight transport can be evaluated by integrating i(x) across 
the surf zone: 



X b 

= / i(x) dx 

* 0 

where x^ is the surf zone width as shown in Fig. (2.3). 

Komar (1977a) calculated the distribution of both the long- 
shore current and the longshore sand transport. Fig. (2.4), 
based on the longshore current analysis by Longuet-Higgins 
(1970) . The distribution shows the maximum longshore current 
at approximately ( x/x ^ = 0.65) and the maximum longshore sand 
transport at approximately (x/x^ = 0.80) . 

Sawaragi and Deguchi (1979) utilized circular bed load 
traps which enable them to measure the on-offshore sand 
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breaker zone 




Figure (2.3). Schematic of defining the location and 

the water depth at the breaker line 
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Figure (2.4). Longshore current and sand transport 

distributions determined by Komar 
(1977a). [from Komar (1977a)] 
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movements as well as the longshore transport distribution 
across the surf zone. The results showed sediment transport 
to be relatively independent of wave steepness H q /L o , but 
dependent on the sand grain size. The maximum sand trans- 
port is shown (Fig. 2.5) to occur at approximately 
(x/x^ = 0.6) and maximum longshore current at ( x/x ^ = 0.4) . 

Bailard (1982) extended his earlier work in which he used 
second order Stokes theory to describe the wave motion. He 
included the Ostenddorf and Madsen (1979) longshore current 
model to predict the distribution of sediment transport rate 
for different ratios of the velocity of the flow to the sedi- 
ment fall velocity (Fig. 2.6). He finds the maximum long- 
shore current at approximately (x/x^ = 0.6) with the maximum 
longshore sand transport at (x/x^ = 0.9) . 

Results from different studies examining the longshore 
current and longshore sand transport distributions are given 
in Table I . 
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X/ X b 



Figure (2.5) . Laboratory measurements of the distribu- 
tion of the longshore current and the sand 
transport rate showing the positions of 
the maximum sand transport and maximum 
longshore current relative to the breaker 
zone. [based on data of Sawaragi and 
Deguchi ( 1979 ) ] 
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DIMENSIONLESS LONGSHORE TRANSPORT RATE. I 







DIMENSIONLESS DISTANCE OFFSHORE, x* 



Figure (2.6). Sediment transport and velocity profiles 

for different values of the velocity ratio 

u mb 

— — , where u ^ is the oscillating water 

velocity magnitude at the breaker line and 
w is the sediment fall velocity. [from 
Bailard, 1982] 
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DIMENSIONLESS LONGSHORE CURRENT. V 



TABLE I 



Location of Maximum Longshore Current (v/v max ) and 
Longshore Sand Transport (q/qmax) Distributions From 

Different Models 



Model and 
Year 


x /*b 


Remarks 


v/v 

max 


q/q 

^max 


Thornton 

(1973) 


— 


1.0 


Field data and 
theoretical relation- 
ship (2.11) with 

B =0.08 
s 


Komar 

(1977a) 


0.65 


0 . 80 


Theoretical rela- 
tionships (2.5) , 
(2.12) 


Sawaragi 

and 

Deguchi 

(1979) 


0.40 


0 . 60 


Laboratory 

measurements 


Bailard 

(1982) 


0 . 60 


0.90 


Theoretical rela- 
tionships 
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III. DESCRIPTION OF LEADBETTER BEACH EXPERIMENT, 
SANTA BARBARA, CALIFORNIA 



A. GENERAL 

Leadbetter Beach, Santa Barbara, California was chosen as 
the second NSTS experimental site, after Torrey Pines Beach, 
California, primarily because it is adjacent to a total sand 
trap forming a spit off the Santa Barbara breakwater, and 
because it has relatively straight and parallel nearshore depth 
contours. The NSTS Santa Barbara experiment began in October, 
1979 with a 14-month sediment trap study measuring sand 
accumulation, coupled with a one-month intensive experiment 
during 27 January through 25 February 1980. 

Exceptional storms were encountered towards the latter 
part of the intensive experiment. During the period 16-20 
February 1980, the weather conditions deteriorated greatly 
approaching the 50-year storm. Significant wave heights in 
excess of 2 meters and longshore current speeds of 150 cm/sec 
were recorded. These conditions coupled with spring tides, 
high winds and wave set-up, eroded the foreshore and back 
beach areas approximately 2 meters in elevation. These 
phenomena coupled with high quality data, make the Santa 
Barbara experiment very unique and significant for near-shore 
transport studies (Gable, 1981) . This study emphasizes the 
very intensive storm period of 17-18 February, 1980. 
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B. COASTAL GEOMORPHOLOGY 



The shoreline between Point Conception and Santa Barbara 
has the unusual (for the California Coast) east-west orienta- 
tion on a predominantly north-south coast. The coastal areas 
of Santa Barbara are composed of marine sedimentary rock and 
Monterey Shale-bluffs 30 m high fronted by sand and cobble 
beaches varying from 0 to 30 m in width (Corps of Engineers, 
1970) . The major sources of sand for the Santa Barbara 
beaches are the Santa Maria and Santa Ynez rivers (see Fig. 
3.1). Some sand comes from the erosion of the coastal cliffs 
(Trask, 1952). Bascom (1951) suggests some contribution by 
wind blown sand, but neither of these two sources are rela- 
tively important. Textural analysis of the sands indicates 
the beaches are composed of fine to medium grain size and are 
relatively well sorted with a mean diameter of 0.236 mm. 

C . WEATHER 

The climate in the Santa Barbara area is classified as a 
Mediterranean type and characterized by warm, dry summers and 
mild, wet winters. The summers are mostly dry because a 
semi-permanent high pressure area covers the Eastern North 
Pacific Ocean, which deflects eastward storms to the north. 
During winters, this Pacific high migrates southward and 
weakens, allowing occasional frontal systems to move through 
Southern California. 

The prevailing winds for the West Coast are from the 
northwest due to the semi-permanent Pacific high pressure and 
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Figure (3.1) . 



Location of 
(USGS, 1976) 



Santa Barbara Channel 
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blow generally parallel to the coastline. However, due to 
the orientation of the coastline at Santa Barbara and the 
protection from Point Conception, the winds are generally 
light and variable from the southerly direction. The aver- 
age wind speed is less than 3 m/sec and the typical pattern 
is that of a breeze. 

Most extratropical winter storms that affect the Santa 
Barbara area move southeastward from the Northeast Pacific; 
on the average, 1 or 2 such storms influence the Santa 
Barbara channel region every year. 

In February 1980, the Eastern Pacific high pressure sys- 
tem was weaker than normal and displaced southeast of its 
normal position. At the same time, the low pressure center 
was also displaced south of its normal position. The surface 
pressures were as much as 18 millibars below normal east of 
the low center. The Eastern Pacific high pressure system 
was almost non-existent, which shifted the storm track suffi- 
ciently to bring storms onshore in Central and Southern Cali- 
fornia (Miller, 1980). 

Two series of frontal systems hit the study site. A 
strong southeasterly occurred from 13 to 16 February and a 
very intensive southwesterly from 17 through 21 February. 
Because these storms were southerly, they produced large 
westerly ocean swells and waves which moved into the Santa 
Barbara channel. The relatively warm water was a major factor 
which intensified the storm and brought much heavier rain 
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than normal. Precipitation was 13.2 cm above normal at this 
time of the year. All these factors, coupled with extreme 
high tides and wave height, caused extreme beach erosion. 

D. WAVES AND WIND CLIMATE 

The wave climate at the experimental site is either locally 
generated inside the Channel Islands or propagated from the 
open ocean. Normally, the predominant waves are swell from 
the west-southwest and the west (240°-270° true) entering 
through a narrow window between San Miguel Island and Point 
Conception with an average wave height of 91 cm and an average 
wave period of 12 seconds (Weigel, 1950). The waves must 
refract almost 90° to reach Leadbetter Beach (Fig. 3.2) . As 
a consequence, the swell waves reaching Santa Barbara are 
greatly reduced in height and arrive at large angles due to 
the severe refraction (O'Brien, 19 50) . 

As a result of the origin of the waves, the offshore 
islands, and the refraction of the waves, the wave pattern 
in Santa Barbara is nearly the same for all times of the year, 
with the exception of the waves from the occasional southeast 
storms. The swell from the open ocean, refracting and break- 
ing at high angle of incidence, create a strong easterly 
longshore current. The occasional southeasterly storms also 
create reversals where the longshore current is to the 
west . 

The local wind waves from the south are usually insignifi- 
cant with only choppy seas and small waves. The local storms 
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Figure (3.2). Refraction diagrams for Santa Barbara and 

vicinity for waves of 10-second period: 

(a) waves from the southwest; (b) waves from 
the south and southeast; (c) expanded refrac- 
tion drawing, waves from the west, west- 
southwest, and southwest; (d) expanded refrac- 
tion drawing, waves from the southeast (after 
O'Brien, 1950) . 
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from the southeast (125°-145° true) have a maximum fetch of 
145 Km and can propagate waves toward the Santa Barbara Coast 
with significant wave height ranging up to 488 cm (Weigel, 
1950). During the storm period (16-20 February, 1980), waves 
generated thousands of kilometers away with unlimited fetch 
were funneled through the wave window of (240°-270° true) . 

The wave period was 14-16 sec with significant wave height 
approaching 200 cm from the southwest and west. The signifi- 
cant wave height was 190 cm from the southeast on February 16 
and 185 cm from the west on February 20, 1980. 

E . EXPERIMENT 

The start of the experiment was timed such that a strong 
negative low tide occurred during the daylight hours allowing 
nearshore instruments to be installed on a dry beach. The 
tides are semidiurnal with a mean range of 110 cm referred 
to mean lower low water. The tides ranged from -10 cm to 
+200 cm (MLLW) during the intensive experiment. 

A plan view of Leadbetter Beach (Fig. 3.3) shows the loca- 
tion of the instruments, the experiment baseline and the five 
range lines. It is noted that the five range lines are not 
equally spaced. 

Series of current meters, pressure sensors and run-up 
meters were used to measure the wave characteristics during 
shoaling, breaking, and run-up. The location of the two wave 
slope arrays is shown in (Fig. 3.4). 
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Figure (3.3). Plan view of the study site showing 

the experimental baseline and the 




Figure (3.4) . 



Location of the two S pressure 

xy 

sensor arrays. J 
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1 . 



Incident Wave Field Measurements 



Two, 6 meter square arrays, composed of four Kulite 
solid state pressure sensors were used to measure an un- 
biased estimate of the spectrum of the radiation stress, S . 

The square arrays are less efficient than a linear array in 
estimating a complete directional spectrum but the narrow 
wave window provided by the offshore islands of Santa Barbara 
Beach produces an acceptable wave climate estimate. The 
location and orientation of the two slope arrays are shown 
in Fig. (3.4). The water depth at both arrays is approximately 
9 m and the sensors were mounted 120 cm above the bottom. 

2 . Profile Measurements 

Profile measurements at Santa Barbara consisted of 
daily wading profiles of the beach face and shallow nearshore, 
coupled with offshore surveys to a depth of 4 m using a trac- 
tor and boat survey for pre- and post-experiment measurements 
to 12 m depth. Five range lines, extending essentially per- 
pendicular to the beach and covering 200 m alongshore, were 
established for the daily wading profile (Fig. 3.3). The off- 
shore distances along range lines are referenced to the distance 
from the experimental baseline (y-axis) . 

The wading profiles were done using conventional rod 
and level survey techniques on a daily basis. Range lines 
were defined by two survey points marked by flags, which a 
rodman aligned before each survey reading. Rod stations were 
measured using a specially constructed plastic-coated steel 
survey line marked at 5 m intervals. An engineer's level and 
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rod provided vertical control. The profile survey was ter- 
minated when the water became too deep for the rodman, or 
the breaking waves made it impossible to plumb the rod. These 
surveys were completed at the lowest possible tides occurring 
during daylight hours. Each survey was closed to the original 
benchmark on land to verify survey accuracy. Data were 
entered on site by hand into a Tektronix computer terminal 
within a short time after each survey. Profiles were plotted 
immediately to visually search for errors which were due to 
either data entry or a bad survey point which could be cor- 
rected by redoing the field survey immediately. 

Generally, small scale changes occurred in the profile 
during most of the intensive month-long experiment (Fig. 3.5) , 
except during the storm where the beach was cut vertically 
2.5 m and horizontally about 60 m within almost a four day period. 
Although it was impossible to maintain the nearshore instru- 
mentation during the storm, the deep water instruments con- 
tinued to work, providing an extremely valuable data set on 
beach erosion. 



40 



FLEVAIiON ABOVE MSl (Ml 




LEADBE T TER beach profiles RANGE no. 3 

28 JAN. TO 25 FEB-. ' 80 



Figure (3.5). Beach profile time history. 
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IV. FIELD DATA ANALYSIS 



A. BEACH PROFILES AND CONSERVATION OF SAND TRANSPORT EQUATION 

Five beach normal ranges at Leadbetter Beach were surveyed 
on a daily basis starting from the established baseline on 
the beach to a distance about 100 m offshore. Standard rod 
and level technique described in (III.E.b) were used. Data 
from the boat-fathometer survey conducted on 25 February, 1980 
for the five ranges were used to extend each profile out to 
the shallow water limit at x = 180 m where (2. < Jj_) . 

Linear interpolation was applied twice, first to calculate 
the missing data points, and second to define five equally 
spaced ranges (50 m apart) . The five ranges were averaged 
for both 17th and 18th of February to produce a daily mean 
profile. Then, a "mean beach" prof ile was defined by averaging 
these two daily profiles. Fig. (4.1) shows the mean beach 
profile, relative to mean sea level (MSL) and indicates no 
well developed bar structures or trough along the profile. 

Based on the difference in elevations between the 17th 
and 18th February surveys at each location, the net longshore 
transport distribution is calculated using the conservation 
analysis of sand transport described later in this section. 

A "mean net longshore transport", <3 mean , distribution was ob- 
tained by averaging the net transport distributions for the 
five ranges. The "mean net sand transport" distribution and 
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Figure (4.1). Mean beach profile gives an average 

profile slope of 0.037 relative to 
HWL, LWL and MSL. 
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Figure (4.2). Sand transport distribution inferred from 

the field data based on the profile changes 
showing a maximum transport at the beach face. 
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the standard deviation calculated at each point are shown in 
Fig. 4.2 indicating a small variance with acceptable scatter 



around the q distribution, 

^mean 



The q distribution goes to zero after a distance of 

^mean 

100 m offshore as shown in Fig. 4.2 because the wading pro- 
files on the 17th and 18th are extended after this distance 
by the same fathometer data taken on 25 February. 

Profile changes between sequential measurements can be 
approximated by a parallelopiped with sides Ax, Ay and the time 
varying depth (Ah) . Considering the equilibrium of the 



parallelopiped element (Fig. 4.3), the net sand trnasport rate 

9q x 

in the x-direction is q Ay = -r-^AxAy and the net sand 

x net d x a 

3q v 

transport rate in the y-direction is q^ net Ax = -r^-AyAx. 



The total net sand transport rate (q T 



Ay + q 



Ax) is 



‘x net J ^y net 
balanced by the rate of the net volume changes over the ele- 



ment giving 



^xnet^ + Snet Ax 



(1-a) ^jr AxAy 



(4.1) 



where a is the sand porosity and is taken as 0.40. Equation 
(4.1) can be written as 



3q x 3q ^ 
3x + 3y 



(1-a) 



3h 
3 1 



(4.2) 



Since a strong longshore current was observed during 17 
and 18 February, (see Fig. (6.3)), a significant longshore 
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Figure (4.3). Conservation of sand transport with 

changes in the beach profile. 
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sand transport was expected. Therefore, it is hypothesized 

that the longshore sand transport is much larger than the on- 

offshore sand transport, i.e., q T7 >> q^. The profile changes 

in Fig. 3.6 do not show any well-defined bar or trough 

structures which suggests the on-offshore sand transport is 

3q x 

negligible (i.e., — = 0). Thus, equation (4.2) becomes 

d X 

Sq y , , . 3h , . _ . 

Ty ~ (1 " a) St (4 ‘ 3) 

Equation (4.3) states simply that the changes in the beach 
profile with respect to time, are proportional to the net 
longshore sand transport. The transport rate/unit width of 
on-offshore distance, q , can be calculated since water depth 
(h) is measured as a function of time (t) . The boundary con- 
ditions are needed to solve (4.3). For the offshore boundary 
condition it is assumed 

q = 0 at h = h 

^y o 

where h Q is the water depth at the shallow water limit, 

h Q = jjj, and L is the wavelength. Given the dominant wave 

period of 14 sec (f = 0.0703 Hz) , h was calculated to be 

p o 

5.0 m at a distance approximately 180 m from the baseline. 

The nearshore boundary condition requires the sand transport 
to go to zero at the shoreline. It appears from Fig. 4.2 
that this condition is satisfied from the data at the baseline 
at x = 0. The net sediment transport rate is calculated by 
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integrating the net cross-shore distribution which is assumed 
to be similar to the longshore sand transport distribution. 

For discussion purposes, two regions within the surf zone 
are defined. A swash region is defined from the baseline to 
the intersection of HWL with the profile, i.e., (from x = 0 
to x - 30 m) ; a seaward region is defined from x - 30 m 
to x - 180 m. The swash region indicates significant 
erosion on the beach face (Fig. 4.2), which will not be considered 
in this study because the model does not include sediment 
transport on the beach face. The seaward region shows the 
eroded profile with a small accreted area around x = 90 m, 
which appears with a negative sign. 

B. RADIATION STRESS, WAVE HEIGHT AND DIRECTION AT OUTER 

BOUNDARY 

Measurements from a 610 cm square, four pressure sensor 

array in 8-m water depth are used to produce an unbiased 

estimate of the spectrum of radiation stress. The pressure 

arrays worked throughout the experiment including the storm 

period. The array measurements are also used to calculate 

the height and the directional characteristics of the waves. 

The radiation stress (S ) is given by (Longuet-Higgins , 

xy 

1962) 



n 

S = / pu w (z,t) v w (z,t) dz (4.4) 

* -h 



where u^ and v^ are the crossshore and alongshore components 
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of wave velocity. The velocities are calculated from the 
pressure gradients measured at the array applying linear wave 
theory spectral transformations. For a particular component 



u j (f ,z,t) 



1 r i cosh k (h + z) , 3p . _ . ~ 

to p cosh k (h + z ) 3x. ' ^ ’ 

m j 



(4.5) 



where the z^ is the measured depth of the pressure sensors, 
a) is the radial frequency and j refers to the two horizontal 
components. The term between the brackets is the spectral 
transfer function H(f) relating the measured pressure gradients 
to inferred velocity components. Substituting equation (4.5) 
into (4.4) , using the measured pressure head ^ and integrating 



S (f) 
xy 



= [- 



P_2 

o? k(cosh^k(h + z ) ) 

m 



Ch „ (f) 

3p 3p 

3x'3y 



(4.6) 



where Cg ^ (f) is the co-spectrum between the measured 

dMy 

pressure gradient components. 

Wave heights were determined from the surface elevation 
records using the "zero-up-crossing method" in which the wave 
height is defined as the difference between the maximum and 
minimum occurring between two consecutive zero-up-crossings. 
Surface elevation records were first linearly detrended to 
exclude the effects of the rising and falling of the tides 
and then high-pass filtered with a cut-off frequency of 0.05 
Hz to exclude surf beat, followed by a low-pass filter with 
a high frequency cut-off (0.3 Hz). 
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The filtering was accomplished by Fourier transforming the 
signals, zeroing the Fourier amplitude coefficients in 
the filtered frequencies, and inverse transforming the 
complex spectrum to obtain the filtered time series. The 
entire 68-minute-record was transformed at one time to mini- 
mize the end effects which result in spectral leakage, and 
to obtain maximum resolution giving very sharp roll-off at 
the filter cut-offs. 

Due to the large number of current meters, current data 
are utilized applying linear theory to infer wave heights 
from the velocities (Guza and Thornton, 1980). Surface ele- 
vation time series were obtained from current measurements 
using the complex Fourier spectra of the horizontal velocity 
component U(f), V(f). The complex surface elevation spectrum, 
X(f), is calculated applying the linear wave transfer function, 
H(f) 



X (f ) = H(f) • V(f) (4.7) 

Then, the complex surface elevation spectrum is inverse 
transformed to obtain the surface elevation time series from 
which the wave height distribution is calculated. Surface 
elevations were also inferred from pressure signals in a simi- 
lar manner by transforming the pressure records using linear 
theory. 

A mean incident wave angle (a) and a representative fre- 
quency (f ) are defined for input into the models such that 
P 



49 



»p 

using a and f to calculate the total radiation stress S , 
P xy' 

the same value is obtained as that measured by the pressure 
sensor array. 

The radiation stress (S ) in equation (4.4) can be ex- 
pressed applying linear wave theory for long-crested waves 



S (f) = E(f) n(f) sin a cos a (4.8) 



The total radiation stress (S ; is calculated using 

xy v 



xy 



= / S (f) df 



xy 



(4.9) 



T T 

Then, the calculated S is equated to its equivalent S 

xy ^ ^ xy 

in (4.8) for monochromatic (f = f ) waves approaching at a 

P 

mean angle a 



xy 



E n(f ) sin a(f ) cos a(f ) 
P P P 



(4.10) 



T 

where E is the total energy. Applying linear theory rela- 
tionships a mean incident wave angle a is defined as: 



a ( f ) 
P 




(4.11) 



where the representative frequency f^ is the peak frequency 
of the energy spectrum. Since the radiation stress is calcu- 
lated using the pressure gradients, the kinetic energy calculated 
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using the pressure gradient is introduced for consistency 
of the analysis. It is assumed the kinetic energy E^ is 
half the total energy. Then equation (4.11) yields 



“<y 



2 sin 



-1 



xy 

E k n( y 



(4.12) 



The kinetic energy is calculated using the pressure gradient 
component spectrum 



E k (f) = |H(f) | 2 [G 3p (f) +G 3o (f)] (4.13) 

9x 3y 

2 

where | H ( f ) | is the linear wave transformation function 
given by 

| H (f)| 2 = PaL sinh(2kh) (4.3.4) 

4kw cosh z k(h + z ) 

m 

The mean incident wave angle (a) at the pressure sensor's 
location was then refracted shoreward to the 5 m contour (at 
the shallow water limit) . A refraction calculation was per- 
formed over a bathymetric chart of the experimental site that 
shows the beach has straight and parallel contours shoreward 
of 5 m depth (Fig. 4.4). Applying linear refraction to the 
measured mean angle (a = 12.4°) at the pressure sensor location 

and peak frequency (f = 0.07 Hz) , waves reach the 5 m contour 

P 

with a refracted angle a = 9.0° where the calculated refrac- 
tion coefficient K equals 0.985. 
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Figure (4.4). Refraction diagram for Santa Barbara 

(17 February, 1980) . 
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The H is calculated using linear refraction theory: 
rms 2 



H 



rms 



H 



rms ( 8m) 



= (■ 



'S 8mi 1/2 , b 8m, 1/2 



( -r> 






1/2 



(K r ) 



(4.15) 



giving a refracted H rms equivalent to .85 m at water depth 

h^ = 5 m, which is used as the initial wave height for the 

model. Refracted wave rays were extended to get the angle 

of approach at the breaker line calculated to give = 6.0‘ 

The breaker line was defined at maximum H . a, was used 

rms b 

to calculate the longshore sand transport using the wave 
power expression. 

Table II gives various calculated parameters of S^, a, 

f , H at different times before, after and during the 
p rms 3 

17th and 18th of February. Orientation of the coordinates 

leads to a negative sign for both S and a, indicating the 

xy 

waves are coming from the west. 
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TABLE II 



Calculated Parameters at West Pressure Sensor 
(During 17-18 February) 



Date 


time 


depth 

(m) 


tide 

(cm) 


S 

XY 2 
(Erg/cm ) 


a ° 


? p 

(HZ) 


H rms 

(cm) 


o 

00 

o> 

H 


0730 


9.2 


+ 71 


-0.43x10 s 


-5.8 


0.094 


65.6 


■8 

fri 

5 


1800 


7.8 


-71 


-1.76x10 s 


-12.4 


0.070 


88 . 8 


r- 

HI 


2230 


9.4 


+ 86 


-2.16xlO S 


-13.0 


.0.062 


96.2 


O 

CO 

o> 

H 


0920 


9.5 


+96 


-2.61xlO S 


-15.9 


0.070 


96.2 


l 

5 

CO 

H 


2200 


9.2 


+ 74 


-1.25xlO S 


-15.0 


0.094 


70.0 
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V. LONGSHORE SEDIMENT TRANSPORT MODEL 



A littoral sand transport model is developed based on 
the energetics approach expressed by Thornton (1973) in which 
the energy transformation is described by the random wave 
model of Thornton and Guza (1983a) coupled with the longshore 
current model of Thornton and Guza (1983b) . The input to 
the predictive model are the offshore rms wave height, mean 
direction and mean period. Both an analytical model for plane 
sloping beaches and a numerical model for arbitrarily bottom 
profile are presented. 

A. DEVELOPMENT OF SAND TRANSPORT MODEL 

Assuming stationary wave conditions with a beach of 

straight and parallel contours, the gradient of the wave 

dECg 

power in the on-offshore direction — - ^ x is balanced by the 
ensemble averaged dissipation due to both wave breaking 
<£j 3 > and bottom friction <£f> (Thornton and Guza, 1983) such 
that 



dECg x 

dx 



<£, > + <£-> 

b f 



(5.1) 



The wave power is expressed using the linear wave theory 
relationships : 



E 



1 

g p g H 



2 

rms 



p g / H 2 P (H) dH (5.2) 
o 
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(5.3) 




£ u + , 2 k h 



2 sinh 2kh 



] cos a 



where H is the root mean square wave height, a is the 
rms 

mean wave direction and k is the wavenumber associated with 



Dissipation due to wave breaking < £j ;) > can be modeled 
after a periodic bore, applying conservation of mass and momen- 
tum at regions of uniform flow up-and-downstream of the bore 
as shown in Fig. (5.1) . The average rate of energy dissipa- 
tion per unit area for single breaking wave is calculated 
(Stoker , 1957) : 



where the wave height H is measured as the maximum to mini- 
mum of the bore, Q is the volume discharge of flow per unit 
area across the bore described by simple linear bore theory 
after Hwang and Divoky (1972) , and B is a breaker coefficient 
of 0 (1) . The coefficient B is an unspecified parameter in 
the model, which is determined from the data. 

The Rayleigh distribution is assumed to describe the 
waves everywhere with the implied assumptions that the waves 
are very narrow band such that they can be described by a 
single frequency f and mean direction a. Thornton and Guza 
(1983) predict the fraction waves which are breaking at 
each location. Using this information, the average rate of energy 



peak frequency f . 



bore 





(5.4) 
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Figure (5.1). Periodic bore used to describe 

spilling breakers. 



(5.5) 



dissipation due to breaking <e^> is given by 



./- B 3 f _ 

3 /it p 

b 16 H y 4, 5 rms 

Y h 



where y is a coefficient relating H to h (taken to be 0.44 

it* 111S 

by Thornton and Guza, 1983) . The average frictional dissi- 
pation <e^> at the bottom boundary layer for the ensemble 
(all waves) calculated using Rayleigh probability distribution 
is : 



<e f > 



P C f 
16/tT 



2tt f H 3 
p rms , 

sinh kh J ' 



(5.6) 



where C^. is the friction coefficient. Thornton and Guza 
(1983a) point out that <e^> is relatively more important for 
gently sloping beaches and for low frequency waves, i.e., the 
shallower the depth the greater the frictional dissipation. 
They neglect the frictional dissipation as it has a small con- 
tribution compared with <£j D > f as will be done here also. 

Substituting equations (5.2), (5.3), (5.5), and (5.6) 

into equation (5.1), yields: 



d_ 1 
dx 8 



p g h 



rms 



C cos a 

g 



3/? 0 9 f p B 

16 4~5 

Y h 



H 



rms 



(5.7) 



H rms can eva l uate< 3 analytically, assuming small inci- 
dent wave angle (cos a - 1) and limiting the analysis to 

ll 1 

shallow water (^- <_ ^-) , which gives 
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(5.8) 



H 



rms 



^ h 9 / 10 [l-h 23 / 4 (- 



23/4 5/2 



) ] 



-1/5 



where 



and 



H 2 h 1 ^ 
o o 



23 ,g,l/2 y tan 8 

a - IF ( ?> b 3 f 

p 



The value of tan 8 is the beach slope and the subscript "o" 
refers to the input conditions at the outer shallow water 
depth limit. When the water depth gets very shallow, (5.8) 
becomes 



H 



rms 



= a 



V5 h 9/10. 



as 



(5.9) 



which indicates that the wave height in the inner surf zone 
is related to the depth and independent of the initial condi- 
tions in deeper water. 

With the assumptions of straight and parallel depth contours 
and stationary wave conditions, the longshore current dis- 
tribution across the surf zone can be estimated using the 
longshore momentum flux balance (Thornton and Guza, 1983) 



3S 



_X£ 

3x 



- T b 



(5.10) 



59 



is the radiation 



where t, is the bottom shear stress and S 

b xy 

stress defined by Longuet-Higgins and Stewart (1962) (for 
linear theory) as 



xy 



_ „ „ sm a 

E C cos a — ~ — 
g C 



From linear wave refraction 



sm a 



sm a 



= constant 



Thus , 



3 S 



xy 

3x 



sm a 



— [i— E C cos a] 
dx g 



(5.11) 



The term in brackets in the R.H.S. is defined in (5.1) 
as <£ b >. Therefore (5.11) reduces to: 



3S 



xv 

— jfc 

3x 



sm a 



<e b > 






(5.12) 



The average alongshore bottom shear stress is defined by 
Longuet-Higgins (1970) , assuming the quadratic law 



t b = P c f |u t | (V + v w )J (5.13) 

where u t is the total instantaneous velocity due to waves 
and current, V is the longshore current, v^ is the longshore 
component of the oscillating wave velocity and c^ is the bed 
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shear stress coefficient. The bar over equation (5.13) indi- 
cates averaging in time over one wave period. Assuming a 
small angle of wave incidence (a < 15°), and a weak current 
case following Liu Ph. L-F and R. Dalrymple (1978) , equation 
(5.13) simplifies to 

x b = p c f |uj V (5.14) 

where u is the mean cross-shore wave velocity, 
w 

|uj is specified using the shallow water wave theory 
relationships and Rayleigh wave height distribution (Thornton 
and Guza, 1983) such that: 




1 

4 



/5ILx 1/2 
1 h ; 



H 



rms 



(5.15) 



Substituting (5.11), (5.12), (5.14), (5.15) into (5.10) 

and solving for the mean longshore current velocity, yields: 



V 



sm a. 



P c 



B 



u 



<£,> 

b 



w 1 



* 1/2 ■ u 6 
f g sin a H 
p o rms 



c f y 



h 



9/2 



(5.16) 



The longshore sand transport distribution is obtained 

by substituting the expressions for the change in energy flux 
3 EC 

( v 9 ) , equation (5.1) and the longshore current equation 

o X 

(5.16) into equation (2.11) to give 
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1/2 



3/2 



(5.17) 



B 



sin a 



g(l | U I 

* p 1 w 1 
s 



P C H 



<£,> 

b 



Substituting (5.7), (5.15) and (5.16) in equation (5.17) 

yields : 



B 



6, 7 



gd -j-) Y b h 
M s 



3 9 19 

sin a f B y H x * 1/2 
^ o _£ rms j 



(5.18) 



Equation (5.18) has the observed characteristics that outside 
the surf zone as h reaches h Q , q s (x) gets small and when h 
goes to zero at the shoreline, q s (x) approaches zero. A 
maximum transport occurs between these two boundaries. 

For a planar sloping beach, the analytical solution for 
H (5.8), which is limited to shallow water, is substituted 

ill O 

into (5.18) to obtain the cross-shore transport distribution. 

For arbitrary bottom profiles (but straight and parallel depth 

contours) , a simple forward differencing scheme is used 

to numerically integrate equation (5.1) for H rms . The 

calculated H value is substituted into (5.18) . The energy 
rms 

flux balance is also numerically integrated from initial 
offshore distance (at depth h Q ) to the shoreline with a 
spatial interval Ax 



E Cg 



x i+1 



= E Cg 



x 1 1 



- Ax<<£k> 



+ <£ f > 



i} 



x 1,2,. «.,n 



(5.19) 
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Equation (5.19) starts at the initial offshore distance 

(i = 1) where H , h , f , a are given. The model calculates 
o o p o 

the forward step by subtracting the total dissipation from 
the calculated energy flux evaluated from the previous step. 
The model predicts H g each step, using linear theory rela- 
tionships (5.2) and (5.3) until i approaches n at the shore- 
line. Then, the predicted H values are used to calculate 

rms 

the longshore current using (5.16). The sand transport dis- 
tribution is evaluated using (5.18). 

Integrating q in (5.18) over limits of offshore distances 
for the study zone, the volume transport rate (Q ) is calcu- 
lated as 



Q z = f q s dx (5.20) 

o 

where x is the offshore distance at which the initial condi- 
o 

tions of waveheight (H ) and water depth (h^) are given. 

The coefficient B in (5.18) is used to calibrate the model. 

s 
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VI. MODEL RESULTS AND COMPARISON WITH PROTOTYPE DATA 



A. MODEL RESULTS FOR H AND V 

rms 

Field data for the Santa Barbara experiment are utilized 

to calculate the parameters S , a , f and H which in 

r xy o p rms 

turn are used as input for both the expressions (5.8) and 

(5.19). Model coefficients y, B and c- were determined from 

H r 

the data. The ratio, y = — was found to be 0.44 from the 

slope of the h vs. H curve. Fig. (6.1) , in the inner surf 
^ rms ^ 

zone. This agrees with the value reported earlier by Thornton 
and Guza (1983) for the Torrey Pines Beach experiment. 

H rms va ^- ues inferred from the pressure and velocity 
measurements at time 1800 on the 17th of February are compared 
with the wave transformation model for various breaker coeffi- 
cients, B. Using an iterative scheme, based on the least 
squares error criterion, B equal to 1.1 gave the best fit to 
H rms ^ ata (Fig. 6.2a) . To check on the generality of the 
model, the waveheight transformation model was run for 5 Feb- 
ruary using B = 1.1. It also gave the best model fit to 
the data (Fig. 6.2b). 

An optimal bed shear stress coefficient, c^ , was calculated 
comparing the longshore current model with the measured long- 
shore current values, V. Since there are only three longshore 
current measurements available on 17 February due to the des- 
truction of instruments by the storm, the data from 5 Feb- 
ruary were again used. With B = 1.1, c^ was varied to 
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rms 



70 




h ( cm ) 



Figure (6.1). Root mean square waveheight vs. water 

depth for 5 February, 1980, giving 
y = 0.44 in the inner surf zone. 
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X (m) 




X (m) 



Figure (6.2). Model H as a function of offshore distance 

compareS m with measurements indicated by x. The 
model is applied to data; (a) on 17 February; 

(b) on 5 February, to obtain B which was 
found to be 1.1 in both days. 
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0.005 



obtain the best fit to the measured V. A value of c f = 
gave the best fit (Fig. 6.3b) . 

B. MODEL RESULTS FOR CROSS-SHORE SAND TRANSPORT DISTRIBUTION 

The empirical coefficient B s in equation (5.18) is used to 
calibrate the model and will be explained later. In all the 
calculations and figures, the offshore distance (x) is 
measured starting from the experimental baseline at x = 0. 

The average net transport between profiles on 17 and 18 
February are to be determined. A weighted time difference 
of 25 hours was found between the two mean profile surveys. 

The model predictions of H , V and q at each x-location 
^ rms' ^s 

depend mainly on the water depth (h) . Therefore, two cases 
are considered to test the sensitivity of the models to the 
variations in the water depth. In the first case, the water 
depth does not vary with time. In the second case, the water 
depth varies by hourly time increments according to tidal 
variation. 

1 . Model Results for Constant Water Depth 

In this case study, the water depths are held constant 
at relative MSL. The numerical model utilizes the actual 
bottom and predicts H , V and q . The model results are 
summarized in Fig. 6.4 indicating the capability of the 
model to predict the increase in waveheight due to shoaling 
and the decrease due to breaking processes. The longshore 
sand transport and the longshore current have a similar dis- 
tribution, which implies that the sand transport may primarily 
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V (cm/sec) V (cm/sec) 



80 




X (m) 




X (m) 



Figure (6.3). Longshore current (68-minute average) as a 

function of offshore distance compared with 
measurements indicated by x. The model is 
applied to data on; (a) 17 Feb., showing only 

a few instruments working; (b) 5 Feb., for which 

= 0.005 was found to give the best fit. 
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be governed by the longshore current (Fig. 6.4). The influ- 
ence of rapid depth changes is reflected in the small pertur- 
bation of the V and q g curves at x = 65 m, which is associated 
with a local flatness of the beach slope, as can be seen in 
the "mean beach profile," Fig. 4.1. The maximum waveheight 
is used to define a "mean breaker line" location, x^, for 
discussion purposes and normalization of the results. The 
analytical model is applied to a plane sloping beach. A slope 
of 0.037 was found as the best fit to the "mean beach" pro- 
file using Fig. 4.1. The normalized distributions resulting 
from the analytical model are summarized in Fig. 6.5, giving 
the maximum sand transport at x/x^ - 0.80, which is similar 
to the results of Komar (1977a) . 

Comparison between the sand transport distributions using 
the numerical and analytical models is given in Fig. 6.6. 

The sand transport peak of the numerical model, which 
uses the actual bathymetry, is shifted offshore compared 
with the analytical distribution using a planar sloping beach. 
It may be seen in equation (5.18) that q^ is proportional 

to (-i,) , which means q increases as h decreases for the 

h 7 3 

same location. The actual bottom profile (Fig. (4.1)) used 
in the numerical model is deeper near the breaker line com- 
pared with the plane sloping beach used in the analytical 
model; therefore, the calculated maximum q g in the numerical 
model is shifted offshore relative to the analytical model. 
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Figure (6.4). Model predictions of H (••••), V ( ) 

and q ( ) for actua5 s bathymetry holding 

depth s constant. 
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Mrmt/maxIHIrm* or V/Vmax or q/qmax 




Figure (6.5). Normalized distribution for H , V and 

q calculated by the analytical model 
(Sonstant depth); H (••••)/ V ( ) 

and q s ( ) . rmS 
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sand transport rate (MU/DAV/M) 

0.0 10.0 30.0 10.0 40.0 10.0 00.0 TO.O 10.0 00.0 100.0 110.0 130.0 




Figure (6.6). Average sand transport distribution 

calculated for constant depth case; for 

analytical model ( ) and numerical 

model ( ) . 
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2 . Model Results for Time-Varying Water Depth 

The models predicting H , V and q were run for 
c rms 

the 25 hours using a one hour time step and correcting the 
depths for the tidal elevations at each time step. The 
resulting 25 sand transport distributions were averaged over 
the 25 hours to give a predicted mean longshore sand trans- 
port distribution between profile surveys on 17 and 18 
February. 

Mean sand transport distributions including tidal 

variation of both the analytical and numerical models are 

given in Fig. 6.7 which shows a good agreement with a peak 

shift towards offshore for the numerical model following the 

previous discussion. When considering the tides, the averaging 

causes the distribution to spread onshore to x = 30 m with 

less q than the first case for the same coefficients. 

^s max 

In both cases the maximum longshore sand transport occurs at 
almost the same location as the longshore current which indi- 
cates again that the sand transport is primarily driven by 
the longshore component in both models. The sediment trans- 
port distribution predicted by the numerical model is not 
smooth, unlike the analytical curve. This is because of the 
use of the actual depth profile. 

C. MODEL RESULTS AND COMPARISON WITH PROTOTYPE FIELD DATA 
The net longshore transport rate distribution for Santa 
Barbara Beach is calculated using the profile changes during 
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•and transport rat* (KflJ/DAV/Kfl) 

0.0 10.0 30.0 10.0 40.0 SO.O 00.0 TO.O 90.0 00.0 100.0 110.0 00.0 




offshore distance (?/;) 



Figure (6.7). Average sand transport distribution _ over 

the daily tidal cycle calculated using 
variable depth; for analytical model 
( ) and numerical model ( ) . 
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17 and 18 February. Most of the estimated sediment trans- 
port occurs at a distance 20 m from the baseline indicating 
a significant beach face erosion, which is beyond the scope 
of the present study. 

Since the volume sand transport rate calculated from the 
data is the net volume change rate in the beach profile, and 
does not represent the total longshore sand transport as the 
model predicts, care must be taken in the interpretation of 
the comparisons. 

A dimensionless sand transport distribution for the field 
data is compared with results from the numerical model in 
Fig. (6.8). The measured maximum net longshore sand trans- 
port occurs on the beach face at x - 20 m. The hypothesis that 
the predicted longshore transport distribution is proportional to 
the net sediment transport distribution, calculated from 
the profile changes does not appear to be the case. Signifi- 
cant profile changes occur at the beach face which are not 
included in the model calculations. It is also possible 
that the sand eroded from the beach face was moved in the 
offshore direction outside 100 m and was not included in the 
measured profiles, which would also invalidate the original 
hypothesis . 

D. THE EMPIRICAL COEFFICIENT B AND LONGSHORE SAND TRANSPORT 

s 

The empirical coefficient B g , as described by Thornton 
(1973) , has a magnitude of order (1) and is equal to 
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q/qmax (DOMENSOONIUESS SAND TRANSPORT) 




offshore distance (.W 



Figure (6.8). Normalized sand transport distribution 

inferred from profile changes vs. 
numerical model results with variable 
depth . 
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( 6 . 1 ) 



B = 7— %T 

s b f h 

where e is the transport efficiency of the available wave 
power and f^ is the dynamic friction factor associated with 
sediment transport. The factor f^ is equal to: 

_ , _ sediment shear stress 

r h ~ normal stress at the bottom ' 



while b is a proportionality coefficient between the ratio of 
sand velocities (^) and the square root of the water particle 
velocity ratio (— ) ^^ . 

Because the data do not provide information about each 
individual parameter, all parameters are combined 
in B s , which is determined empirically by model calibration. 

The total longshore sediment transport rate (Q^) in 
(m'Vday) is calculated by integrating the area under the sand 
transport distribution 



< 2 * 



x=180 

/ q s dx 



The calculation based on profile changes does not give the 
total longshore sand transport but gives the net transport 
which is much less in this case than the longshore transport 
predicted by the model. Therefore, equation (2.5) suggested 
by CERC (1977) is used as a reference to calibrate the long- 
shore sand transport model. 
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The calculations of the longshore sand transport rate 
require the wave condition at the breaker line H ^ , h^ 
and a^. Field data give H rins ^ equal to 0.9 8 m. Fig. (6.2a) 
and h^ to be 2.96 m knowing the beach slope besides the angle 
c*k equals six degrees from the refraction diagram. Fig. (4.4) . 

3 

calculated using equation (2.6) gives 4600 m /day, which 
requires B g to be 0.18. 

A summary of the results if given in Table III. The 
difference between the two cases of depth (constant and 
varying) for the analytical model, was found to be insignifi- 
cant. A similar result was found for the numerical model. 

The calculated using the numerical model is less than 
calculated from the analytical model, because the numerical 
model used the actual bottom profile. The numerical and 
the analytical models gave essentially the same maximum 

cross-shore transport distribution q for the same case 

^ ^max 

of either constant or varying depth, Fig. (6.6) and Fig. 

(6.7) . 



Johnson (1953) estimated the longshore sand transport 
at Santa Barbara Beach to be 220,000 m'Vyear, or about 600 
m'Vday on the average. The Q, predicted by the model gives 
almost eight times the above estimated average daily rate. 

The high value of predicted by the model reflects the 
storm conditions, which caused significant longshore transport. 
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VII. SUMMARY AND CONCLUSIONS 



The principle objective of this study is to predict 
the longshore sediment transport distribution. Special 
emphasis was placed on testing an analytical and numerical 
model against field data acquired from an experiment at 
Leadbetter Beach, Santa Barbara during the intensive storm 
period 17 to 18 February, 1980. 

Field data show a strong longshore current coupled with 
significant changes in the beach profile with no evidence 
of on-offshore sand movement (bars or troughs) along the 
profiles. Therefore, it is hypothesized that the changes 
in the profiles are due to the variations of sand transported 
alongshore by the observed longshore current, and that the 
on-offshore sand transport is negligible. The net sand 
transport calculated from the profile changes is expected 
to have a distribution similar to the cross-shore sand 
transport distribution. 

Initial values for the models of H , a and f are 

rms o p 

o 

obtained from the field data. Model coefficients y = 0.44, 

B = 1.1 and c^ = 0.005 are determined by model fitting. The 
developed model is capable of predicting the increase in wave 
height due to shoaling and the decrease after wave breaking. 
The model predicts the longshore current (V) and the longshore 
sand transport (q g ) at each location. A maximum q g is 
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predicted at distance x = 0 . 8 x^ which agrees with Komar 
(1977a) . The experimental transport coefficient, B , is found 
to be 0.18 instead of 0.08 by Thornton (1973). The maximum 
longshore current, V max , occurs at approximately the same 
location as (q ) , indicating that the sand transport 

distribution is primarily determined by the longshore current 
distrib ution. 

The net change in the beach profiles used to infer the 
net sand transport distribution, shows q s maximum at the 
beach face. Since the sand transport model does not consider 
the swash region, the net profile change distribution does 
not agree well with the predicted cross-shore distribution, 
and the original hypothesis does not appear to be valid. As 
a result, sand transport mechanics in the swash zone appear 
to be important in determining why the maximum sand trans- 
port occurs at the beach face. The swash is driven primarily 
by infragravity waves (long waves) due to surf beat and edge 
waves which have not been included in the surf zone models 
where only the depth-limited sea-swell frequency band of 
waves are considered. Therefore, further studies need to 
model the swash zone in order to describe the onshore water 
wave boundaries in a more realistic manner. 
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